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1, INTRODUCTION 

Most of electronic devices experienced thermo-mechanical fatigue during services. Soldered joints 
are used to connect different electronic parts and therefore provided the way to transfer electrical currents. 
Due to the imperfection of joints, cracks may appeared especially at the interface between two materials. 
In general, the concept of Linear Elastic Fracture Mechanics (LEFM) is used to characterize such cracks in 
order to determine the reliability of the joints. For example Huang et al. [1] provided a general methodology 
for calculating mixed mode stress intensity factor of interfacial cracks. In this work, only a single crack is 
focused considering single-edge crack in plain strain condition. Tsai [2] performed a study on the mode I 
energy release rate of soldered joints containing a single crack using finite element analysis considering a 
method of virtual crack closure. It is only considered a single crack between two parts. Another similar work 
can be found in [3]. Fundamental work on the interfacial cracks between two dissimilar materials are 
available in [4-7]. Most of these papers discussed on the effect of dissimilar materials on the single crack. 
It is hard to find the works on the multiple cracks or crack interactions in different material behavior. 

Comninou [8] provided an overview work on the interfacial cracks between two different materials 
considering the LEFM. Atkinsons’ first and second models are applied in analyzing the interfacial cracks. 
In this work, only a single crack is considered. There are several works found in discussing the role of 
multiple interfacial cracks for an example Monfared et al. [9] studied on the dissimilar materials with 
multiple cracks. It is positioned exactly between two materials where the crack geometry is identical. 
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On the other hand, Bidokhti and Shahani [10] also studied the multiple cracks however the cracks 
are non-aligned. Meanwhile, Wen et al. [11] also studied the multiple surface cracks under fatigue loading. 
In this work, cracks are co-planar with both identical and non-identical sizes of three-dimensional models. 
The results are then compared with ASME, BS7910, API579 and GB/T19624 codes of standard. Then, they 
suggested a new combination rule and concluded that it always yields a reasonable estimation with necessary 
conservatism for various initial crack depth, material constant and relative sizes of two cracks. Ismail [12] 
investigated numerically the multiple crack interaction in bi-materials plate under mode I tension. Two plates 
of different materials are perfectly bonded assuming that there is no chemical reaction at the interface. 
Two cracks are positioned at a distance of b for the interfacial boundary. In this work, identical crack length 
is used to study the crack interactions when different mechanical mismatch is used. The effect of mechanical 
mismatch can also be found in [13, 14]. Azlan and Ismail [13] studied the inclined crack in bi-material plates 
and Ismail [14] examined the crack eccentricity at different distance measured from the interfacial line. 
However, both works focused on the effect of mechanical mismatch on the single crack. Based on the 
literature found above, it is seemed that the study on the crack interaction when the crack sizes are un- 
identical are hard to find. 

Therefore, this paper numerically presents the role of multiple cracks embedded at the joint. 
There are two cracks of un-identical crack length and embedded at the joint material. ANSYS finite element 
program is applied through the use of ANSYS Parametric Design Language (APDL). Two important 
parameters are used to study the crack interactions such mechanical properties mismatch and crack length. 
J-integral is used to calculate the crack driving forces and then convert into stress intensity factors. 
It is selected since the use of singular element around the crack tip is unnecessary. 


2. RESEARCH METHOD 

In this work, three multiple cracks are modelled using ANSYS APDL finite element program as in 
Figure 1. It 1s assumed that the model fulfilled a plain strain condition and subjected to axial tension. 
Since the model is symmetrical, only a quarter model is considered and modelled as in Figure 2 where the 
vertical left side is assigned as a symmetrical line and on the other hand the bottom line is also assigned as a 
symmetrical side except the crack faces. 

Figure 3 shows the quarter finite element model of multiple crack under tension force. It is meshed 
using PLANE183 element type assuming to fulfill the plain strain condition. The model is constructed using 
ANSYS Parametric Design Language (APDL). In this work, J-integral is used as a crack driving force before 
it is converted into Stress Intensity Factor (SIF) as in (1). Since J-integral is used, it is not required to use 
singular elements around the crack tip. Figure 4 shows the deformed plate contained multiple cracks under 
tensile force. In order to study the effect of dissimilar materials, the factor of mechanical mismatch, 
e 1S used as in (2). Three values of e are used such as 1, 5 and 10. In this work, the height ratio of joint is kept 
constant where h/W = 0.02 and Table 1 summarized the crack parameters used. 





D (1) 


where, J is the value of J-integral, K is a stress intensity factor, E is a modulus of elasticity and v is a 
Poisson’s ratio. 


ee od (2) 


where, e is a mismatch factor, E; is a modulus of elasticity for plate material and E> is a modulus of elasticity 
for soldered material. Before any further analysis, it is important to validate the present model with the 
existing solutions. It is hard to find the multiple cracks in dissimilar or solder joints. Then, it is assumed that 
the solution obtained from homogenous plate and single crack (Lan et al. [7] and Yan & Miao [15]) as in 
Figures 5 and 6. It is observed that for both type of cracks (double edge and central cracks), the present 
cracks are well agreed with the existing solutions and verified for the further applications. 

Stress intensity factors, K obtained from ANSYS finite element program is normalized as follow: 
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where, F is a normalized SIF or geometrical correction factors, K is the SIF, ois the applied stress and a is 
the crack length. In order to study the interaction among the cracks, an interaction factor, vis then introduced 
as below: 


F 
aa (4) 


edge crack 


where Feentrat_crack Ad Fredge crack are the normalized SIFs obtained from central crack (crack length, a;) and 
edge crack (crack length, a2), respectively. 








H 
Figure 1. Multiple three crack at the solder joint Figure 2. Quarter model of the solder joint with 
multiple cracks 
Figure 3. Finite element model of multiple cracks Figure 4. Finite element results of multiple cracks 


Table 1. Summary of Crack Parameters 


a,|/W a>/W WwW h/W a) Qa? 
0.1 0.1 50 1 5 5 
0.2 0.1 50 I 10 i 
0.3 0.1 50 I 15 = 
0.4 0.1 50 I 20 i 
0.5 0.1 50 I 25 a 
0.6 0.1 50 I 30 fs, 
0.7 0.1 50 I 35 5 
0.8 0.1 50 I 40 5) 


*Note: The parameter a,/W 1s also varied starting from 0.1 to 0.8. 
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3. RESULTS AND ANALYSIS 

This section is divided into two separate sections. The first section deals with the effect of 
normalized edge crack length, a./W on the interaction factor, v The second part on the effect of mechanical 
mismatch, e on the interaction factor, v. 


3.1. Effect of Normalized Edge Crack Length, a2/W on the Interaction Factor,y. 

Figure 7 shows the roles of normalized central crack length, a)/W on the interaction factor, y when 
mechanical mismatch is kept as a constant parameter. As shown in (4), the factor of crack interaction is 
defined as the ratio between SIFs of central cracks over the edge cracks. Higher the interaction factor 
indicating that the SIFs of central cracks are higher than the SIFs for edge cracks. It is observed that as in 
Figure 7, interaction factors increased as the function of the length of central cracks. Many research works 
have reported similar results and can be found in [7-14]. 
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Figure 7. Effect of normalized edge crack length, a)/W on the SIFs when different mechanical mismatches, e 
are varied, (a) e = 1.0, (b) e = 5.0 and (c) e = 10.0 
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On the other when the normalized edge cracks increased the interaction factors are also decreased. 
This is due to the fact that for the case of E; > E>, the resistance of joint is lower than the upper and bottom 
plates. This condition created a large joint deformation and therefore producing higher values of SIFs. 
In comparison these SIFs are higher than the SIFs when E; < E>. Figures 7(a)-7(c) reveal the interaction 
factors for different values of mechanical mismatch, e = 1.0, 5.0 and 10.0, respectively. It is observed that the 
mechanical mismatch played an important role in determining the interaction factors especially for the case 
of a,/W = 0.1 and 0.2 where increasing the mismatch values capable to reduce the interaction factors. 


3.2. Effect of Mechanical Mismatch, e on the Interaction Factor, +». 

Figure 8 shows the effect of mechanical mismatches, e on the interaction factor, y when normalized 
edge cracks, a,/W are varied. In general, when the value of mechanical mismatch increased, there is 
insignificant effect on the interaction factors especially when higher value of a2/W is used. On the other hand, 
lower interaction factors can be observed when longer edge cracks are introduced. This is indicated that the 
interaction factor is significantly related to the relative length between two interacted cracks. Figure 6(a) 
shows the behavior of interaction factor for the case of a./W = 0.1 where both cracks are similar in size. 
When both cracks are identical, the effect of mechanical mismatch on the crack interaction is obvious where 
it is decreased when the mismatch factor is decreased. 
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Figure 8. Effect of mechanical mismatches, e on the interaction factor, y for different normalized edge crack, 
(a) a2o/W = 0.1, (b) ao/W = 0.2 and (c) ao/W = 0.3 


4. CONCLUSION 
Based on the numerical simulations, several main conclusions can be listed as below: 
a) It is observed in general that the interaction factor reduced when the mechanical mismatch is increased. 
b) It is also observed that the interaction factor converged into a single value when the a2/W is greater 
than 0.3. 
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